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The field-theoretical model describing multicritical phenomena with two coupled order parameters with n\\ and 
n± components and of 0(n\\)<sO{n±) symmetry is considered. Conditions for realization of different types of 
multicritical behaviour are studied within the field-theoretical renormalization group approach. Surfaces sep- 
arating stability regions for certain types of multicritical behaviour in parametric space of order parameter 
dimensions and space dimension d are calculated using the two-loop renormalization group functions. Series 
for the order parameter marginal dimensions that control the crossover between different universality classes 
are extracted up to the fourth order in e = 4-d and to the fifth order in a pseudo-e parameter using the 
known high-order perturbative expansions for isotropic and cubic models. Special attention is paid to a particu- 
lar case of 0(1)® 0(2) symmetric model relevant for description of anisotropic antiferromagnets in an external 
magnetic field. 
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1. Introduction 

The concept of universality plays a paradigmatic role in the modern statistical physics. Accordingly, 
continuous phase transitions can be grouped into universality classes (see, e. g. . Systems within the 
same universality class are characterized by the same set of critical exponents governing the scaling 
behaviour of their thermodynamical functions. Therefore, one of the aims of a theoretical description of 
a system is to establish its universality class. 

In the theory of critical phenomena it is standard now to use methods of field theoretical renormaliza- 
tion group (RG) . Within these methods, a stable fixed point (FP) corresponds to the universality class. 
For systems with complex internal symmetries described by <p 4 theories with several couplings, several 
different nontrivial FPs may exist. Depending on global parameters of a system, these FPs can interchange 
their stability causing the system to trigger from one universality class to another. The lack of a stable 
FP can even mean that a continuous phase order transition is transformed into a discontinuous. These 
global parameters (that effect the FP stability) are spatial dimension d and the dimension n of the order 
parameter (OP). In the rc-d-space, the regions of stable FPs are separated by borders and the n{d) curves 
define the OP marginal dimensions that control the crossover between different universality classes. 

In this paper we are interested in the stability borders and marginal dimensions for a model with 
two coupled OP fields, namely, the model with O(ny) ® 0(n±) symmetry . Such a model describes, 
amongst other systems @], anisotropic antiferromagnets in an external magnetic field jloHl^l . 

Conditions for realization of different types of multicritical behaviour, that are defined by the relation 
between the dimensions of the OPs n«, n±, were obtained already in the first nontrivial approximation 



jf s nu, n±, we: 

of the field-theoretical RG for d < 4 Oj [lj, [l8J] . They determine the stability regions in the parametric 
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n\\—n± plane for three FPs: isotropic HeisenbergYV of 0{n\\ + n±) symmetry, decoupled FP at which OPs are 
ordering separately, and biconical FP. The two-loop studies in d — 3 show qualitatively similar results 
Hill , although significantly changing the quantitative picture in n\\ - n± plane. Five-loop results for three- 
dimensional 0{n\\) ® 0(n±) model (20fl confirm the obtained picture, producing only slight corrections. 

Since the previous studies of multicritical behaviour in the 0[n\\) ® 0{n±) system concentrated on 
d — 3 case, in this paper we consider the dependence of marginal dimensions of 0{n\\) & 0(n±) model on 
space dimension d. Our motivation is caused by the fact that even a small change in d can produce crucial 
effects on the critical behaviour, in particular, changing the universality class of a system. The rest of the 
paper is organized as follows: In section[2]we present the 0{n\\) ® 0{n±) model and its RG description. 
Then, our aim is to analyse the conditions for realizing different scenarios of multicritical behaviour. In 
section [3] we present the results obtained within the two-loop approximation based on the e-expansion 
as well as on the fixed d approach. We devote the next section|4]to the results for marginal dimensions of 
0{n\\) ® 0(?2j_) in higher order approximations. We end the paper with section [5] where our conclusions 
are presented. 



2. The model and RG picture of its multicritical phenomena 

The model with 0{n\\) ® 0(n±) symmetry can be obtained from the well-known 0(«)-symmetrical 
model 12111 . splitting its n-component OP (pQ into two: cp±o and cp\\o that act in orthogonal subspaces with 
dimensions n\\ and n±, respectively (n\\ + n± = ri): 

0o = K ±0 )- (2.1) 



lio 

Then, separating the Ginsburg-Landau-Wilson functional of 0(n) symmetry one can present the effective 
Hamiltonian of the 0{n\\) © 0{n±) model in the form: 

d 1 i d 



d d x\- r_L0_Lo • 4>io + - £ V^io • V,-(/>_lo + - hHo • <P\\o + - £ ^iHo • Vi<p\\o 
L z 1 i=i z 1 i=i 

+ -^j 1 (0±o-0±o) +^(0||o-0no) +^(0±o-0±o)(0i!o-0i!o) 



(2.2) 



where three couplings u±_ and & x should be introduced instead of the only one in the 0{n) symmetric 
model, and f± and f\\ are connected with the temperature distance to the critical line for cf>±o and 0no, 
correspondingly. 

The first mean-field analysis of the model with two coupled OPs was performed in order to describe 
the supersolids fiill (see also El). It shows that the character of the multicritical point in such a phase 
diagram depends on the sign of u± u\\ - u 2 x . For a positive sign, a tetracritical point is realized, while for 
a negative sign, it is a bicritical point. Going beyond the mean field theory, fluctuations should be taken 
into account. This is achieved by the field-theoretical RG approach | 2], in which the large-scale behaviour 
of the system is connected with the stable FP of the RG transformations. The transformation of the fourth 
order couplings {11} in <2.2> under renormalization is described by /3-functions. 

The /3-functions for 0[n\\)@0{n±) model were known in a one-loop approximation [8]. The next order 
approximation has been found in the massive flill as well as in the minimal subtraction RG schemes 1 16] . 
In the minimal subtraction scheme, the ^-functions were also calculated in the five-loop approximation 
Eoll . although explicit expressions were presented only for 0(3) ® 0(2) symmetry j23ll . Here, we work 
with /3-functions obtained in two-loop order flill within the minimal subtraction RG scheme jj^Eill : 

(n±+8) 9 n» , (3rc_i_ + 14) o 5nii , n« q 

Pu L = -eu_i + — -u 2 ,+—u 2 x -—-± -u\ -u L ul-—ul, (2.3) 

6 6 12 36 9 

(«l+2) («||+2) 2 2 (nx + ray+16) 3 

p Ux = —EU X H U±U X H UxU\\ + -Ux Ux 

6 6 3 72 

(n x +2) 2 («|! +2) 2 5(w ± +2) 2 5(«n+2) 2 

uzui utu\\ u 1 Ux u x u„, (2.4) 

6 6 " 72 1 72 11 

(My +8) , n± , (3wij+14) o 5n± , n± , 

Put = -EU11 + U7< H UZ UX, UllUx ut ■ (2.5) 

y 11 11 6 11 6 12 11 36 11 9 
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Here, {u±,u x ,u\\} - {u} are renormalized couplings and the space dimension d enters the /3-functions via 
e = 4- d. 

The FPs {u*} of the RG transformation are found from the zeros of the /3-functions 

Pu,[{u}) = Q (2.6) 

with i =_L, ||, x. A stable FP possesses positive eigenvalues a>i, (1)2, 0J3 (or their real parts) of stability 
matrix dfiilduj. 

The stable FPs for 0{n\\) ® 0(ra_J are already known from the one-loop studies For d < 4 and for 
sufficiently low OP dimensions satisfying 

n± + n\\ < 4, (2.7) 

only the isotropic Heisenberg FP J€ of 0(n± + n\\) symmetry with {u^ — u* = up is stable. When n± (or 
n\\) increases breaking H2.7t . still with 

n±n\\ +2{n± + n\\) < 32, (2.8) 

FP 36 interchanges its stability with biconical FP £8 {u* ± + m* + ut}. For values of n± and n\\ that are 
above the condition <Z8l , FP SS looses its stability, while the decoupled FP @ {u* ^ 0, m* = 0, wjf ?i 0} be- 
comes stable. According to these one-loop results, the multicritical behaviour of the 0(1) © 0(2) model 
is governed by FP J€ (connected with bicriticality) for all space dimensions d < 4. However, within the 
higher order calculations, the stability of FPs depends not only on ray , n± but also on d. Using resum- 
mation procedures for the two-loop RG functions at d = 3, one can show that the conditions of the FPs 
stability <2.7> and (278) are drastically shifted to smaller values of OP components Si . In particular, 
in the case n\\ = 1, n± = 2 FP M (connected with tetracriticality) appears to be stable in a two loop order 
. Resummation of higher orders c-expansion [ 20] does not change this result. 



3. Stability border-surfaces within a two-loop order approximation 

As noted above, the stability of FPs @, 88, J€ is dependent on three parameters n\, n± and d. There- 
fore, the borders between regions for which one or another FP is stable, form surfaces in the parametric 
space n\\ - n±- d: f[n\\,n_i,d) = 0. We call them border-surfaces (BSs). 

Two alternative ways are used in practice to analyze RG functions and to get universal quantities, 
in particular, marginal dimensions. In one approach, i.e., the e-expansion, the solutions are obtained 
as a series in e and then they are evaluated at the value of interest (for instance, at e = 1 for d = 3 
theories). Alternatively, one may fix the space dimension d to a certain value and directly solve a system 
of non-linear equations obtaining the FP coordinates numerically (2^1 . In the next two subsections we 
use these approaches to obtain marginal dimensions of the 0{n\\) ® 0(n±) model within a two-loop RG 
approximation. 

3.1. BSs from e-expansion 

We start our analysis with establishing the border between the stability regions of the decoupled FP 
& and the biconical FP S&. As it was noted in I16I1 . two of the FP @ stability exponents correspond to 
the stability exponent of the O(n) model w^'"': wf = w jei - n t\ tof = w^'" 1 ', while the remaining one is 
defined by 

wf = dp Ux lu x \g. (3.1) 

Since w 1 ^'"' is always positive, only wf governs the stability of the FP @, changing its sign depending 
on ray, raj_, d. Therefore, the surface between stability regions of FPs ® and 38 can be extracted from the 
condition of J3.lt vanishing. Substituting the £-expansion for the FP Qs coordinates into J3.lt one collects 
terms up to e 2 and sets the result equal to zero: 

[(13ran +44) (ray +2) (ra ± +2)(13ra ± +44) 
£ [ 2(ra,| +8) 3 2(ra ± +8) 3 



"II 



n± -4 



2(ran+8) 2(ra_i+8) 



(3.2) 
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(a) (b) 



Figure 1. (Color online) BSs between different universality classes of the 0(n\\) e 0[n±) model obtained 
(a) by applying an e-expansion and (b) by using a resummation procedure at a fixed d to two-loop RG 
functions. The left hand (lower) surface separates the stability region of FP J€ (on the left from the sur- 
face) and FP 8& (on the right from the surface). The right hand (upper) surface separates the stability 
regions of FP SM (on the left from surface) and @ (on the right from surface). The vertical line shows the 
position of a system with n\ = 1, n± = 2. The disc on the line indicates the position at d = 3. 



This is analytically solved for e = e{n\\,n±). The result is shown as the right hand surface in figure[T](a). 

The BS between the regions of stability of the FPs SB and J€ can be derived from the condition that 
FP J€ changes its stability. Only one of the three eigenvalues of the stability matrix changes its sign in the 
region considered. Calculating this eigenvalue up to the e 2 order we get the equation for the surface: 

[-5(«l -fray +8) 2 + 66(rc_L + ray +8)-360le I 12 \ 

J + 1=0. (3.3) 



(n± + n\\ +8^ \ nj_ + n\\ + 8 

The surface is also shown in figure[T](a) (the lower left hand surface). 

The limiting borderlines in the plane e — (d — 4), are equivalent to the case when the one loop order 
inequalities <2.7> . <2.8> are transformed into equalities, from which one obtains 

m 2(16- n±) a. 
»f(»x) = r— » nf («_l) = -«_l + 4. (3.4) 

The vertical line in figure [T] (a) presents a system with n\\ — 2, n± — 1, indicating which FP governs 
the multicritical behavior of this system with the change of e. Note that the FP SB is stable in the region 
0.51 < e < 1.04. We are interested in this case, since it describes anisotropic ferro- and antiferromagnets 
in space dimension d = 3. 

3.2. BSs from resummed /^-functions 

Another way to obtain the BSs, is to calculate them from the /3-functions J2.3M2.5> fixing d at cer- 
tain values. Since the RG expansions have divergent |2j nature, the special resummation techniques are 
needed to get convergent results 12911 . The two-loop /3-functions <2.3) - f2\Sl /3 = /3 Ui have a form of polyno- 
mials in renormalized couplings: 

j8({u})= £ c ljk u} x u\ul. (3.5) 
We first represent <3.5> in the form of a resolvent series lioll in one auxiliary variable t : 

F{{u},t)= £ c ijk u[uiu k x t i+j+k ~ 1 = E a a {{u},{c})t a , (3.6) 



where the expansion coefficients a a in B.6t explicitly depend on the couplings and on the coefficients 
Cijk <3.5> . Obviously, F{{u], 1) = fi{{u}). We resume the function H3.61 as a single variable function using 
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the Pade-Borel technique [ 31] and writing its Borel image as: 

a n t a 



F B {t)= £ 



(3.7) 



Analytical continuation of the function (377) is achieved by representing it in the form of a Pade approxi- 
mant fiill . In our case, we use the diagonal Pade approximant [1/1]: 



F a (f)^[l/l](r). 

Finally, the resummed function is obtained via an inverse Borel transform: 



F res = J [l/l](f) e - f . 



(3.8) 



(3.9) 



Applying the above procedure to the two-loop /3-functions J2.3H2.5t at a fixed d and searching for 
their FP solutions with u* = and together with expression d3.1> , where wf = 0, we find a BS, separating 
the stability region of the FP 3> from the FP S& stability region. It is depicted as the upper surface in 
figure [l](b). Searching for the FP solutions with u* L — u* = u,T at which the determinant of the stability 
matrix vanishes we derive the BS between the stability regions of the FPs J€ and §8. This is the upper left 
hand surface in ngure[T](b). 

It is technically difficult to extract the data from the resummed function in the limit e — ► 0. Therefore, 
we present BSs for 0.002 =£ e =£ 1.2, and n\, n± in the range from -0.56 to 5. Limiting borderlines in 
the plane £ — described by <3.4> give us the one-loop (thin) borderlines of figure 1 of , while the 
intersections of the surfaces with the plane £ = 1 give the two-loop (thick) borderlines of figure 1 of Hill. 




Figure 2. (Color online) Locations of the FPs 38 (triangles), J6 (discs) and ® (squares) for n± = 2, ny = 1 
and e changing from 1 (end right hand marks) to with a step size of 0.1. Arrows show the direction 
in which e decreases. Results are obtained using resummation 13.6H3.9) . The line shows the track of a 
stable FP. Dotted part of the line indicates that FP 38 is stable, while the solid part indicates that FP J€ is 
stable. Black disc at the origin indicates the Gaussian FP. 

Similarly to what we did it in figure [l] (a), we present in figure [T](b) the line indicating the stability 
regions for FPs of the 0(1) ® 0(2) model. In this approximation for the /3-functions, FP S& is stable in the 
region 0.66 < e < 1.06, in particular at £ — 1 {d — 3). Let us check how the FP picture changes along the 
line with an increasing d. Varying the space dimension d from 3 to 4 with stepsize 0.1 we can observe the 
drift of FPs 3$, J€, @ towards to the Gaussian FP. The tracks are shown in figure [2]by triangles, discs and 
squares indicating the change of locations of FPs J€ and Q> in u± - u\\ - u x space with the change of d. 
Numerical values of the coordinates of these FPs are listed in table[l] FP S& is stable up to the intersection 
of traces of FP 88 and FP J€, which happens at d « 3.34, where it interchanges its stability with FP X . 
Thus, only the FP J€ is stable starting from the intersection point and up to the Gaussian FP. This would 
mean that in higher space dimensions and at n\\ — 1, n± = 2, the phase diagram contains a bicritical point 
instead of a tetracritical point. 
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Table 1. Coordinates of the FPs SM, 3€, and @ at n\\ = 1, n± = 2 depending on e as found from the resumed 
two loop /^-functions. 



£ 


uf* 


"II 


uf'* 




uf* 


uf-* 
"II 


-1 

1, 


1 1 9 r 7 r 7 
1.1Z / / 


1 9Q r 7/l 
l.Zo / *± 


U.oUlo 


1.UU1D 


1,111 J 




0.9 


0.9112 


1.0039 


0.5273 


0.8434 


0.9569 


1.0971 


0.8 


0.7313 


0.7739 


0.5799 


0.7026 


0.7939 


0.9063 


0.7 


0.5834 


0.5934 


0.5518 


0.5771 


0.6496 


0.7386 


0.6 


0.4596 


0.4502 


0.4863 


0.4650 


0.5214 


0.5906 


0.5 


0.3543 


0.3348 


0.4048 


0.3647 


0.4075 


0.4600 


0.4 


0.2634 


0.2405 


0.3183 


0.2749 


0.3061 


0.3444 


0.3 


0.1843 


0.1627 


0.2323 


0.1944 


0.2158 


0.2420 


0.2 


0.1149 


0.0982 


0.1497 


0.1223 


0.1354 


0.1513 


0.1 


0.0539 


0.0446 


0.0720 


0.0578 


0.0637 


0.0710 



4. High loop order results for marginal dimensions 

The marginal dimensions of the 0{n\\) ® 0{rijJ models can be defined based on the high order RG 
results for simpler isotropic and cubic models. In particular, exact scaling arguments fish connect the FP 
® stability with the critical exponents of the 0(«_J and 0{n\\) models: 



(!)■) — 

2 2 



a{n±) a{n\\) 

+ — 

v{n L ) v(«||) 



1 1 

(4.1) 



v(rej vC«||) 



where a{n) and v[n) are the heat capacity and correlation length critical exponents of the 0{n) model. 

As it was indicated in jioll . the stability of the FP J€ is defined by the marginal dimension n c of 
the cubic model. Since in the FP J€, the RG functions depend only on the combination n = n± + n\\, the 
resulting marginal dimension can be presented in the form nf{n\\,E) = n c {e) - n\\ . 

In the following two subsections we present an analysis of the marginal dimensions rf[n\\,e), 
nf{n\\,E) based on the five-loop minimal subtraction series for the RG functions of isotropic 13411 and 
cubic models lisll . as well as for the case d — 3 based on the six-loop series for these models lia I37I1 
obtained within the massive scheme 1 38, 391. 



4.1. Five-loop ^-expansions for marginal dimensions 

Let us start with the calculation of nfn\\,e). Substituting the five-loop e-expansions of the 0{ri) theory 
fiill into d4.lt and putting u>f equal to zero, we get the equation for the BS. Keeping n\\ as a parameter 
and expanding in e, we get (n\\, e) in the following form: 

nxin\\,eHn$ +2) = 2(16-n||)-48£+8[3C(3)(« 2 +34?in+10o] + [n 2 +58n l |+148]]i? 2 i|| e 2 

-920aiJ -528n 2 +21376n N + 51584] /3 
- 4 (l06592 + 64480rcn + 13548rc 2 + 1258nj| + Yln\ j C(3) /3 
+ [l8(l00 + 34«n +n 2 ]c(4)-40[550 + 163?i|| +7n 2 j<'(5)/3] J R^ | 2 |i?* || e 3 

+ (3nJj i + 170^-43120«]|-442864^-2069072?i 2 - 4512896«u -3457280] /6 
-50(550+ 163nn +7n 2 ]((6)i?^/3 
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-K(3)(l816192+5011904rc||+3131936rc 2 +623376n 3 +41740?2 4 + 14! 



«|-»h)/ s 



4 (l51424 + 131552«|| + 15728rc 2 - 1250ra 3 - 17«J - £ 2 (3)i?~| /3 
(l06592 + 64480«|| + 13548rc 2 + 1258^ + 17 wj] ((4)R~f 
(4822640 + 2331088nn + 416334rc 2 + 51745rc 3 + 1103n|] <T(5)i?~ 2 /9 



520 + 31792nn + 5826rc 2 + 535rc 3 + 17rc 4 ) £ (7)i?~ 2 /2 



where _R„ = (n + 8) . Expressions for certain physical values of tt\\ are less cumbersome: 



n ± (l,e) = 10. - 16.£ + 22.84224e 2 
«x(2,e) = 7. - 12.£ + 17.76523e 2 - 
ra±(3,e) = 5.2 - 9.6e + 14.43837e 2 



-44.06758c 3 + 113.6428c 4 , 
-34.18402c 3 + 84.07657c 4 , 
- 28.00490c 3 + 67.23923c 4 . 



(4.2) 



(4.3) 
(4.4) 
(4.5) 



The obtained e-expansion diverges, as it can be seen from the growth of the expansion coefficients in 
||4~3)-{|4~5) as well as it follows from the Pade table (H for n±{l, 1): 



Bx(l,D = 



10.0000 
-6.0000 
16.8422 
-27.2253 
86.4175 



3.8462 
3.4092 
1.7981 
4.5288 
o 



3.4773 
3.9728 
2.8846 

o 

o 



2.4576 
3.1128 

o 

o 

o 



9.6637 A 
o 







(4.6) 



The element MAT of the table 14.6) is the value of n±[l, 1) given by the [M/N] Pade approximant at e = 1. 
Here and below, symbol o denotes the approximants which can not be constructed within the order of 
perturbation theory considered here. Usually, the best convergence of the results is observed along the 
main diagonal and the closest sub-diagonals of the Pade table |32j] . However, it appears that the value of 
W_l(1, 1) given by the Pade-aproximant [2/2] differs from those given by [1/2] and [2/1] by an order of one, 
leading to an uncertainty of the numerical estimate. 

To obtain a reliable estimate of wj_(l,e) we rely on the Pade-Borel resummation described in sub- 
section \32\ We obtain a resolvent series by a substitution e — » et. For the obtained expression we build 
the Borel-image, then approximate it by the [3/1] Pade approximant. Performing an integration of the 
inverse Borel transform, we arrive at the result shown in figure[3]with a solid line. In a similar way, we 



4.0, 
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2.0 



























—1 










S 








7 


<— 






4 


► 





1-5 

0.0 0.2 0.4 0.6 0. 



1.0 1.2 1.4 



Figure 3. (Color online) The dependence of the marginal dimensions n® (1) (solid line) and «^°(1) (dashed 
line) on e. The results are obtained based on the five-loop expansion for isotropic and cubic models using 
Pade-Borel resummation with [3/1] Pade approximant (see the text). The diamond denotes the location of 
the three-dimensional O(l) ® 0(2) system. 



get 
[35] 



et n^{n\\,e) using the available five-loop c-expansion for the marginal dimension of the cubic model 



4-2e- 



5 

12 



+ 503^ 2 + (15({4) 25C(5) 



5<T(3) 



2 ) [8 
(15({4) 125C(6) 11515f(7) 



3 8 72 J 

3155£(5) 229(T(3) 2 93f(3) 



I 32 



12 



384 



1728 



144 



128 



1 

3"84 



-U 4 . (4.7) 
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The dependence of nf^{\,e) = n c {e) - 1 on e is obtained as above by the Pade-Borel resummation with 
[3/1] Pade approximant. The result is shown in figure[3]with a dashed line. 

As it can be seen from figure [3] that FP for the 0(1) ® 0(2) model is stable in the region 0.84 < e < 
1.36. The value of n± - 2, d - 3 (denoted by a diamond) is located very close to the boundary n±(l,e). 
Thus, one concludes a very slow approach to the FP. Measurements in 0(1) ® 0(2) systems may show an 
effective critical behaviour with the exponents close to the 0(3) case (l6fl. Anyway, in recent Monte-Carlo 
simulations of the Heisenberg ferromagnet with uniaxial exchange anisotropy, only a bicritical point with 
Heisenberg symmetry was obtained 140(1 . 



4.2. Marginal dimension for of d = 3 in a six-loop order 



Fixing the spatial dimension to d = 3, we can analyze (nj) using pseudo-e expansions (for details 
see I33jl). Introducing the pseudo-e parameter t into 6-loop RG functions of the massive scheme at a fixed 
d — 3 13711 for the 0{ri) model, one can derive critical exponents in the form of pseudo-e expansions and 
substitute them into equation J4.H . Similar to the former subsection, we extract the pseudo-e expansion 
for w® (W||) and present as an example 



"±(D 
n±{2) 

BJ.C3) 



10 - 10.66667T + 5.13069^ 



■2.30752t 3 + 1.69527t 4 



1.98282T 3 



7-8T + 3.99473-T 
5.2-6.4t + 3.24817t 2 



1.76429t 3 + 1.140396t 4 
1.46248T 3 +0.84832T 4 



1.20818T 3 , 
0.84314T 5 . 



(4.8) 
(4.9) 
(4.10) 



The pseudo-e expansions have better convergent properties, as it is known from other studies jisllillliill . 
This is also seen from the coefficients in the series <4.8) - (4~T0t , as well as from a comparison of the Pade 
table <4.6> with the one that follows from the pseudo-e expansion d4.8t : 



m(i) 



10. 


4.8387 


3.7156 


3.2882 


3.1541 


3.0391 


-0.6667 


2.7977 


2.8683 


3.0805 


2.0645 





4.4640 


2.8724 


2.7743 


2.9854 








2.1565 


3.1338 


3.0073 











3.8518 


2.9379 














1.8690 


















(4.11) 



However, the convergence of the results might be spoiled if a pole in the denominator of a Pade 
approximant appears. We demonstrate this below by the Pade table for «j_(3): 



n±(3) = 



5.2 


2.3310 


1.6662 


1.3945 


1.2670 


1.1823 


-1.2 


0.9546 


1.0319 


1.1042 


0.9528 





2.0482 


1.0397 


-0.4474 


1.0623 








0.5857 


1.1226 


1.0705 











1.4340 


1.0112 














0.5909 


















(4.12) 



where by small digits we indicate a result for the approximant [2/2] with a pole for t = 0.944. 

From the Pade-Borel procedure with [4/1] approximant, we get: «®(1) = 2.981. An estimate for nj^(l) 
readily follows from the known result obtained based on the six-loop pseudo-e-expansion n c — 2.862 1 36] . 
Subtraction of 1 leads to the following result n^{\) - 1.862. 



5. Conclusion 

In the present paper we have studied the conditions under which different types of multicritical be- 
haviour are realized for the 0(?i||) © 0{n±) model. These types are related to the three FPs (J€, Q>, 3S), 
and their stability defines the regions in the space of the dimensions of the OPs as well as in the spatial 
dimension where the corresponding multicritical behavior manifests itself. Using the e-expansion for the 
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two-loop /3-functions obtained in the minimal subtraction scheme we derived the BSs separating these re- 
gions. We obtained similar BSs applying the resummation procedure. In the particular case of 0(1) ® 0(2) 
symmetry, we confirm the previous studies finding that the biconical FP associated with a tetracritical be- 
haviour is stable for the case d = 3. In higher space dimensions, the 0{n\\ + n±) symmetrical FP associated 
with the bicritical behaviour is stable. 

Our analysis also made use of the results of higher order approximations within the field-theoretical 
RG approach. At this stage, there were used the scaling arguments connecting the stability of the FPs 
of 0(H||) ffi 0(rcj_) model with the universal quantities of the 0[ri) and the cubic models. Exploiting five- 
loop expressions for the 0{ri) model, we derived an £-expansion for the marginal dimension n®{n\\,e) 
separating the regions of stability for the FPs @ and SS. Applying the resummation procedure to this 
result, we have analyzed the dependence of n®(l,e) on e. Exploiting the five-loop expressions for the 
cubic model we obtained the value of n^fl.e) separating the regions of stability for the FPs J€ and 
S3. Finally, we complete our results by three-dimensional estimates of «®(1) and n^{Y) based on the 
pseudo-£ expansions derived within a six-loop RG approximation. 

These results are also important for the critical dynamics I43U47I1 . The type of a dynamical FP in such 
systems depend, of course, on the static FP values. In order to extend our results to the dynamics of 
antiferromagnets in an external field, further work is necessary. One has to extend this analysis to the 
statics of the corresponding model C l4^ - [5oll . 

Acknowledgement 

This work was supported in part by the FP7 EU IRSES project N269139 "Dynamics and Cooperative 
Phenomena in Complex Physical and Biological Media". 

References 

1. Domb C, The Critical Point: A Historical Introduction to the Modern Theory of Critical Phenomena. Taylor & 
Francis, 1996. 

2. Brezin E., Le Guillou J.C., Zinn-Justin J., In: Phase Transitions and Critical Phenomena, eds. Domb C, Green M.S. 
Academic Press, London, vol. 6, 1976, p. 125. 

3. Amit D.J., Field Theory, the Renormalization Group, and Critical Phenomena. World Scientific, Singapore, 1989. 

4. Zinn-Justin J., Quantum Field Theory and Critical Phenomena. Oxford University Press, 1996. 

5. Kleinert H., Schulte-Frohlinde V., Critical Properties of ip 4 -Theories. World Scientific, Singapore, 2001. 

6. Nelson D., Kosterlitz J.M., Fisher M.J., Phys. Rev. Lett., 1974, 33, 813; doi 10.1103/PhysRevLett.33.813 

7. Aharony A., Bruce A.D. , Phys. Rev. Lett., 1974, 33, 427; doi 10.1103/PhysRevLett.33.427 

8. Kosterlitz J.M., Nelson D., Fisher M.J., Phys. Rev. B, 1976, 13, 412; doi|loTll03/Phys RevB.13T412"l 

9. Aharony A., J. Stat. Phys., 2003, 110, No. 3-6, 659; doi |10.1023/A:1022103717585l 

10. Shapira Y., In: Multicritical Phenomena, NATO ASI Series B, eds. Pynn R., Skjeltrop A., Plenum Press, New York 
and London, vol. 106, 1983, p. 35. 

11. Rohrer H., Phys. Rev. Lett., 1975, 34, 1638; doi |10.1103 /PhysRevLett.34.1638] 

12. Rohrer H., Gerber Ch., Phys. Rev. Lett., 1977, 38, 909; doi 10.1103/PhysRevLett.38.909 

13. King A.R., Rohrer H., Phys. Rev. B, 1979, 19, 5864; doi |10.110 3/PhysRevB.19.5864l 

14. Butera R. A. et al., Phys. Rev. B, 1981, 24, 1244; doi 10.1103/PhysRevB.24.1244] 

15. Ohgushi K., Ueda Y, Phys. Rev. Lett., 2005, 95, 217202; doi |10.110 3/PhysRev Lett.95.217202| 

16. Folk R., Holovatch Yu., Moser G., Phys. Rev. E, 2008, 78, 041124; doi |10.1103 /PhysRevE. 78.041124| 

17. Lyuksyutov I.F., Pokrovskii V.L.., Khmelnitski D.E., Sov. Phys. JETP, 1975, 42(5), 923-926. 

18. Aharony A., In: Phase transitions and Critical phenomena, eds. Domb C, Green M.S., Academic Press, London, 
vol. 6, 1976, p. 357-424. 

19. Prudnikov V. V., Prudnikov P. V., A. A. Fedorenko, JETP Lett., 1998, 68, 950; doi |10.1134/1.567959| [Pis'ma Zh. Eksp. 
Teor. Fiz., 1998, 68, 900 (in Russian)]. 

20. Calabrese P., Pelissetto A., Vicari E., Phys. Rev. B, 2003, 67, 054505; doi 10.1103/PhysRevB.67.054505 

21. Brezin E., Le Gulliou J., Zinn-Justin J., Phys. Rev. B, 1974, 10, 892; doi 10.1103/PhysRevB.10.892 

22. Liu L., Fisher M.J., J. Low Temp. Phys., 1973, 10, 655; doi |10.1007/BF00655458| 

23. Hasenbusch M., Pelissetto A., Vicari E., Phys. Rev. B, 2005, 72, 014532; doi[lb.ll03/PhysRevB.72.014532 

24. 't Hooft G., Veltman M., Nucl. Phys. B, 1972, 44, 189; doi |10.1016/0550-3213(72)90279-9| 



43001-9 



M. Dudka etal. 



25. 't Hooft G., Nucl. Phys. B, 1973, 61, 455; doi 10.1016/0550-3213(73)90376-3 In the minimal subtraction scheme 
only pole terms appearing in vertex functions are subtracted. It is used together with dimensional regularization. 

26. Dohm V., Z. Phys. B, 1985, 60, 61; doi |10.1007/BF01312644| 

27. Schloms R., Dohm V., Europhys. Lett., 1987, 3, 413; doi |10.120 9/0295-5075/3/4/005 

28. Schloms R., Dohm V., Nucl. Phys. B, 1989, 328, 639; doi 10.1016/0550-3213(89)902233 

29. Holovatch Yu., Blavats'ka V., Dudka M., von Ferber C, Folk R., Yavors'kii T., Int. J. Mod. Phys. B, 2002, 16, 
4027; doi:10.1142/S0217979202014760 (A review on application of resummation technique to divergent RG ex- 
pansions.) 

30. Watson P.J.S., J. Phys. A: Gen. Phys., 1974, 7, L167; doi |10.1088/0305-4470/7/18/001| 

31. Baker G.A., Jr., Nickel B. G., Meiron D.I., Phys. Rev. B, 1978, 17, 1365; doi |ToTl 03/PhysR evBT7.1365| 

32. Baker G.A., Jr., Graves-Morris P., Pade Approximants, Addison-Wesley: Reading, MA, 1981. 

33. Le Guillou J.C., Zinn-Justin J., Phys. Rev. B, 1980, 21, 3976; doi |10.1103/PhysRevB.21.3976| 

34. Kleinert H., Neu J., Schulte-Frohlinde V., Chetyrkin K.G., Larin S.A., Phys. Lett. B, 1991, 272, 39; 
doi |10.1016/0370-2693(91)91009-K| Erratum: Phys. Lett. B, 1993, 319, 545. 

35. Kleinert H., Schulte-Frohlinde V., Phys. Lett. B, 1995, 342, 284; doi 10.1016/0370-2693(94)01377-0 

36. Folk R., Holovatch Yu., Yavorskii T., Phys. Rev. B, 2000, 62, 12195; doi 10.1103/PhysRevB.62. 12195 

37. Antonenko S.A., Sokolov A.I., Phys. Rev. E, 1995, 51, 1894; doi [miT 03/PhysRevE.51.1894 

38. Parisi G., In: Proceedings of the Cargrese summer School, 1973 (unpublished). 

39. Parisi G., J. Stat. Phys. 1980, 23, 49; doi 10.1007/BF01014429 

40. Selke W., Phys. Rev. B, 2011, 83, 042102; doi |10.1103/PhysRevE.83.042102| 

41. Dudka M., Holovatch Yu., Yavorskii T., J. Phys. Stud., 2001, 5, No. 3/4, 233. 

42. Holovatch Yu., Ivaneyko D., Delamotte B., J. Phys. A, 2004, 37, 3569; doi |TaT0 88/0305-44 7073"7/ll/002| 

43. Folk R., Holovatch Yu., Moser G., Phys.Rev.E, 2008, 78, 041125; doi 10.1103/PhysRevE.78.041125 

44. Folk R., Holovatch Yu., Moser G., Phys.Rev.E, 2009, 79, 031109; doi |10.1103/PhysRevE.79.031109 

45. Folk R., Holovatch Yu., Moser G., J. Phys. Stud., 2009, 13, 4003. 

46. Folk R., Holovatch Yu., Moser G., Euro. Phys. Lett, 2010, 91 46002; doi |10.1209/0295-5075/9l /46002 

47. Folk R., Holovatch Yu., Moser G., Phys.Rev.E, 2012, 85, 021143; doi 10.1103/PhysRevE.85.021143 

48. Halperin B.I., Hohenberg P.C., Ma Sh.-K., Phys. Rev. B, 1974, 10, 139; doi|l0Tll03/PhysR evBlai39l 

49. Brezin E., De Dominicis C, Phys. Rev. B, 1975, 12, 4954; doi 10.1103/PhysRevB.12.4954 

50. Folk R., Moser G., Phys. Rev. Lett., 2003, 91, 030601; doiil0.1103/PhysRe vLett.91.0 30601 



TpaHI/IMHi BMMipHOCTi A/151 MySlbTI/IKpi/ITI/IMHI/IX <|>a30BI/1X 

nepexofliB 

M. AyAKfP, P. cpo/ibKP, K). ro/iOEiaJ 11 , I". MosepP 

iHCTHTyT (|)i3HKW KOHfleHCOBaHWX CHCTeM HAH YKpaiHH, By/1. CBEHL4il4bKOrO, 1, 7901 1 ZlbBiB, YKpaiHa 

iHCTHTyT TeopeTH4HOi c(>i3HKH Yi-iiBepcuTeTy MoraHa Ken/iepa, A-4040 J1 i H14, ABCTpin 
iHCTHTyr <j>i3HKH Ta 6ioc|)i3HKW YHiBepcuTeiy, A-5020 3a/ibLj6ypr, ABCrpm 

Po3r/iflflaeTbC« TeopeTHKO-no/ibOBa MOfle/ib, 140 onncye My/ibTUKpnTUHHi HBuma i MaeflBa 3B'fl3aHi napawieTpw 
nopaflKy 3 «|| i n± KOMnoHeHTaMM Ta 0{n\\)®0{n±) ci/iMeTpiio. Y paMKax TeopeTHKO-no/ibOBOi peHopMrpy- 

nn BMBHaK)TbCfl yMOBH pea/li3aL4ii pi3HWX TUniB My/1bTWKpHTW4H0i nOBefliHKH. BHK0pHCT0ByKJ4M flBoneT/ieBi 

peHopMrpynoBi djyHKL4iY, po3paxoByiOTbCfl noBepxHi, 1140 po3fli/iflfOTb 06/iacTi CTMKOCTi fl/ia neBHHX TwniB Kpn- 
tuhhoT noBefliHKM, b napaMeTpn4H0My npocTopi BUMipHOCTeti napaMeipiB nopaflKy Ta npocropoBOi BHMip- 
HOCTi d. BuKopucTOByio™ po3K/iafln pj\ft i30TponHOi Ta Ky6i4HOi MOfle/iefi, BiflOMi y bucokhx nop^flKax Teopii 
36ypeHb, OTpuMyiOTbCfl pflflu pj\n rpaHH4Hoi BMMipHOCTi napaMeTpa nopaflKy, flKa KOHTpo/iioe KpocoBep Mix 
pi3HHMH K/iacaMn yHiBepca/ibHOcri, flo 4eTBepToro nopflflKy 3a e = 4-d Ta p,o n'HToro nopaflKy 3a nceBflo-e na- 
paMeTpoM. Oco6/iMBa yBara npufli/iaeTbca BunaflKy 0(l)ffiO(2) cnMeTpn4Hoi MO,qe/ii, flKa B/iacTUBa fl/ifl onucy 
aHi30TponHHx aHTUcfiepoMarHeTHKiB y 30BHii±iHbOMy MarHiTHOMy no/ii. 

KniOMOBi c/ioBa: My/ibTUKpuni^Hi HBi/ima, rpaH^'-iHi BUMipHOcri, peHopMrpyna 
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